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Abstract 

A class of exact steady and unsteady solutions of the Navier Stokes equations in cylindri- 
cal polar coordinates is given. The flows correspond to the motion induced by an infinite disc 
rotating with constant angular velocity about the z-axis in a fluid occupying a semi-infinite 
region which, at large distances from the disc, has velocity field proportional to (x,—y, 0) 
with respect to a Cartesian coordinate system. It is shown that when the rate of rotation 
is large Karman’s exact solution for a disc rotating in an otherwise motionless fluid is re- 
covered. In the limit of zero rotation rate a particular form of Howarth’s exact solution 
for three-dimensional stagnation point flow is obtained. The unsteady form of the partial 
differential system describing this class of flow may be generalized to time-periodic equilib- 
rium flows. In addition the unsteady equations axe shown to describe a strongly nonlinear 
instability of Karman’s rotating disc flow. It is shown that sufficiently large perturbations 
lead to a finite time breakdown of that flow whilst smaller disturbances decay to zero. If the 
stagnation point flow at infinity is Sufficiently strong the steady basic states become linearly 
unstable. In fact there is then a continuous spectrum of unstable eigenvalues of the stability 
equations but, if the initial value problem is considered, it is found that, at large values of 
time, the continuous spectrum leads to a velocity field growing exponentially in time with 
an amplitude decaying algebraically in time. 
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Contract No. NAS1-18605 while the first author was in residence at the Institute for Computer Applications 
in Science and Engineering (ICASE), NASA Langley Research Center, Hampton, VA 23665. Additional 
support was provided SERC. 
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1. Introduction. 


Our concern is with the strongly nonlinear instability of the boundary layer on a rotating disc, the 
important feature of our investigation is that the stability problem which we discuss corresponds to an exact 
solution of the Navier Stokes equations. In addition our investigation uncovers a class of exact steady and 
unsteady Navier Stokes solutions relevant to the flow over a rotating disc immersed in a three-dimensional 
stagnation point flowfleld. 

The importance of the stability problem for the boundary layer on a rotating disc is due to the fact that 
the flow can be thought of as a prototype stability problem for boundary layer flows over swept wings. Some 
years ago Gregory, Stuart and Walker (1955) showed that the flow over a rotating disc is highly unstable 
to an inviscid ’crossflow’ instability associated with the highly inflectional velocity profiles which occur in 
certain directions. Their calculations pointed to the particular importance of a stationary mode of instability 
associated with an effective velocity profile having an inflection point at a position of zero flow velocity. The 
structure of the Gregory, Stuart and Walker mode in the nonlinear regime was later discussed by Bassom 
and Gajjar (1988), the latter authors show that in that regime a nonlinear critical layer structure develops. 
This mode is apparently the one most preferred in an experimental situation if the background disturbance 
level is sufficiently small. However, though most experimental investigations of the rotating disc problem 
have clearly identified the stationary crossflow structure described by Gregory, Stuart and Walker (1955), 
some experiments have pointed to the existence of a second type of vortex structure associated with some 
type of subcritical response caused by another type of instability, see for example Federov et al (1976). 

A possible cause for this second type of stationary vortex structure is the viscous stationary crossflow 
vortex identified numerically by Malik (1986) and described using essentially triple-deck theory by Hall 
(1986). Later MacKerrell (1987) was able to show that this mode is destabilized by nonlinear effects and 
therefore might cause the subcritical instability observed experimentally. However a key feature of the 
mechanism described by Hall (1986) is that the crucial balance of forces leading to instability is one between 
Coriolis and viscous forces, thus in swept wing flows this mechanism is possibly not operational. An alter- 
native source of finite amplitude instability in more general three-dimensional boundary layers is the one 
discussed in this paper. Though we shall formulate and solve the resulting nonlinear interaction equations 
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in cylindrical polar coordinates it is easy to see the relevance of the structure we find to flows more naturally 
described in Cartesian coordinates. 

In Section 2 we shall formulate the nonlinear interaction equations describing the flow over a rotating 
disc immersed in a three-dimensional stagnation point flow. In Section 3 we shall discuss some steady 
equilibrium states of these equations, in particular we describe the nonunique nature of the solutions of 
these equations. In Section 4 we concentrate on the unsteady forms of the interaction equations and discuss 
the linear and nonlinear instability of the flow over a rotating disc in an otherwise still fluid. Our calculations 
in that Section point clearly to a threshold amplitude response of the flowj thus a sufficiently large initial 
disturbance causes an unbounded velocity field to develop after a finite time. The singularity structure 
associated with this ’blow-up’ is discussed in Section 5. Since the interaction equations which we derive in 
Section 3 are obtained without neglecting any terms in the Navier Stokes equations the singularity discussed 
in Section 5 is a singularity of the full Navier Stokes equations in three dimensions. Finally in Section 6 we 
draw some conclusions. 

2« The equations for combined disc-stagnation point flows. 

With respect to cylindrical polar coordinates (r, 6 } z) the Navier Stokes equations may be written 

(~t\ _1 

— + (u • V)u + as = — Vp+J/Au +v ( -jSf-f ^ , 

V 0 ) p \ 0 J (2.1) 

div y = 0, 


where (u, u, tu) is the velocity field corresponding to (r, 0, z) and p, p and u are the fluid pressure, density and 
kinematic viscosity respectively. The operators V and A appearing in (2.1) are the gradient and Laplacian 
operators in cylindrical polar coordinates. We define dimensionless time and axial variables T and £ by 


T=Kt, C = 



(2.2a, b) 


2 



where K is a constant. We seek a solution of ( 2 . 1 ) in the form 

u = Km(C,T) + \{K^U^T)e ii9 + c.c.}, 

v = Orfl(C, T) + \iiKjUiC, T)e ™ + c.c}, 

(2.3o, b, c, d) 

w = (Ki/)$V>((,T), 

?■ = -\K 2 r 2 + vKp{(,T) + ±-y 2 K 2 r 3 {J{T)e™ + c.c.}. 

p 2 2 

Here 7 , A are constants whilst c.c. denotes ‘complex conjugate’. We note that (2.3) reduces to Karman’s 
solution for the flow over a rotating disc if we choose Q = K and set 7 = 0. If we substitute (2.3) into 
( 2 . 1 ) we find the crucial result that the nonlinear terms in the radial and azimuthal momentum equations 
generate no terms proportional to e i4t *. This means that (2.3) is an exact Navier Stokes solution and we 
find that the equations to determine the functions appearing in (2.3) are 

Or + O 3 + 7 J |^| 2 - Q0 = -A + 

Dy -f* 2 HD -f* tDD^ = j 

iDt + 

2tH- tZ)( = 0, 

and 

U T + 2 VU + ViU( - - J + U u . ( 2 . 5 ) 

From now on we shall restrict our attention to the situation when fl = K\ note that Karman’s solution 
is retrieved by letting 7 -» 0. Before writing down boundary conditions appropriate to (2.4), (2.5) it is 
perhaps worthwhile commenting on the motivation for the choice of the special form (2.1). In a recent 
paper, Balakumar, Hall and Malik (1990), investigated the instability of Karman’s solution to nonparallel 
Rayleigh modes of wavenumber n in the azimuthal direction. These high Reynolds number modes can be 
made nonlinear in the manner suggested by the vortex-wave interaction structure of Hall and Smith (1990). 
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In that structure the amplitude of the nonparallel mode with azimuthal wavenumber n is adjusted until it 
drives a mean flow correction comparable with the unperturbed state. For 0(1) values of n it turns out that, 
because of the comparable size of the three disturbance velocity components in the analysis of Balakumar, 
Hall and Malik (1990) , only the n = ±2 modes can be made strongly nonlinear in the manner described by 
Hall and Smith (1990) and then the appropriate form of the disturbed flow is (2.1) . However the structure 
(2.1), suggested by the interaction described by Hall and Smith, is applicable at all Reynolds numbers rather 
than just at high Reynolds numbers where the work of Balakumar, Hall and Malik (1990) applies. Thus we 
can interpret (2.1) as a ‘mean field* type of disturbed flow with the 0-dependent part representing a wave 
superimposed on Karman’s solution which now evolves in time as the disturbance develops. We note that 
in this evolution procedure no terms in the Navier Stokes equations have been neglected. 

We close this section with a discussion of the boundary conditions required to enable a solution of the 


system (2,4), (2.5). We assume that as the flow evolves the mean (with respect to 6) part of the velocity 
field, i.e. (fi, 9, iD) , satisfies 

0=0, 0=1, ID = 0, C = 0, 


0=0, 9 = 0, ( — ► oo. 

Next we assume that the wavelike part of the flow satisfies 


(2.6a, b ) 


u = 0, C = 0, U = e iNT , C~*oo 


(2.7a, b) 


where N is a constant dimensionless frequency. In order that the O, and U equations are consistent with the 
above conditions A and J must be choosen such that 


A = - 7 3 , J = -iNe iNT . (2.8a, b) 

Having made the above choice of boundary conditions we can seek solutions of (2.4), (2.5) which are periodic 
2ir 

in time with period — , the steady states of (2.4), (2.5) are then found by setting N = 0. It follows from the 
form of the nonlinear terms in (2.4), (2.5) that periodic solutions have (ft, 9, t&) independent of time so that 
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(o, V, tB) and V satisfy 

iNU + 2017 + mU ( = U (< + iN 

O 3 + 7 3 |17| 3 - P 3 + tPO< = 7 3 + Q« 

200 + OiP ( = P« 

► 

20 + = 0 

0 = 0, 0=1, tl = 0, 17 = 0, c = 0 

0=0 = 0, < — » oo, 17=1, < = oo J 

where we have replaced £f(£,T) by 17(£) . The nonperiodic solutions satisfy 

Uj> + 2o{7 + HiU^ — U(( + %N e’^ 3 " 

Or + o 3 + 7 3 [C7| 3 - 0 3 + tP0< = 7* + 0«. 

0, + 200 + lP0< = 0«, 

20 + = 0 

► 

0=0, 0=1, to = 0, U- 0, C = 0 
O = 0 = 0, c = 0, U = e iNT , ( = OO 


(2.9) 


( 2 . 10 ) 


0=6(0, 0 = 0(C), U = 17(0, T = 0 J 

Thus the periodic solutions can be found by integrating an ordinary differential system, (2.9), wheras the 
unsteady modes satisfy a parabolic partial differential system. For that reason we have been required in 
(2.10) to give initial conditions to completely specify the problem for 0,0, iD and U . Furthermore we note 
that (2.10) can be regarded as the appropriate nonlinear initial value instability problem for the periodic 
problem (2.9). In the next section we shall discuss the solutions of (2.9), the solution of (2.10) will be 
discussed in §4. 
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3 Equilibrium solutions of the interaction equations. 

In order to begin the solution of (2.9) by some appropriate numerical method it is convenient to discuss 
a limiting form of that system which can then be used to begin the calculation. The first limit we consider 
is 7 -+ 0 in which case the problem for fX t V and i 2 ) becomes uncoupled from that for U and is in fact simply 
Karman’s solution. Thus we know that in the limit 7 ^ 0 , ti'( 0 ) ~ .50, 0'(0) ~ -.61. Another known flow is 
found in the limit 7 — ► 00 . In that limit we write 


n = yu 0 + ..., 

9 = % + ..., 
tD = 7*®o + •••, 

U=U 0 + 


(3-la, b, c, d) 


where the functions appearing in the above expansions are functions of the stretched variable ( = In 
this limit the problem for tlo, Uq, Dio decouples from Do and we obtain the system 

iNU 0 + 2 VqUo + wqUq = CT" + iN, 

+ \U 0 \ 3 + *0^ = 1 + Off, 


2 flo + ti)' 0 = 0 , 


> 


(3.2) 


Do = 0, Dio = 0, U 0 = 0, { = 0, 

fl<) = 0, Uq = 1, £ — > OO. j 

In the special case N = 0 we can relate (3.2) to a special case of the three-dimensional stagnation point flows 


considered by Howarth (1951), Davey (1961), Banks and Zaturska (1988). We recall that, with respect to 
Cartesian coordinates x, y, z Howarth identified the class of exact Navier-Stokes solutions given by 

u = V = Ij-g'iv), w = - (^) V/W + Vg(v)), (3.3) 


where l is a length and 



(3.4) 


Here /', g' satisfy 


r -(/+«»)/" =i+r 

g ,2 -(g + -f)g" = i + -g‘ 

a a 


(3.5) 
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with OL — VjV and subject to 


f = g = f' = 9' = o, r] = o 


(3.6) 


f = g’ = 1, T|-too. 

These solutions correspond to a three-dimensional stagnation point flow above the plane z — 0. In the special 
case a = —1 we can relate /' and g* to the functions Uq, Uq with N = 0 by writing 


2iio = /' - -V>o = / - g , 217 0 = /' 4- 


Davey (1961) gives /" = 1.2729, y" = —.8112 which suggests that for large 7 the solutions of (2.9) with 
N = 0 are such that 


tt'(0) = 1.0427* 


U’( 0) = .23l7*+---. 


(3.7a, 6) 


However we shall see below that the solution of (2.9) is not unique so Davey’s solution correspond to only 
one of our solutions at large values of 7. For nonzero values of jV a similar asymptotic structure can be 
obtained but the coefficients in the expansions (3.7a,b) will, of course, be functions of the frequency. Before 
giving the results of our numerical investigation of (2.9) we note that, for 7 >> 1, the dominant terms in 
the steady state solution of (2.3a,b,c) are such that 


(u, r,tu) ~ rK"y{HQ + TJqcqs29 } —Uq sin 20, 



(3.8) 


where without any loss of generality we have taken JTo to be real. If we transform (3.8) to Cartesian 
coordinates we obtain a velocity field 


u ~ K 7[xGo + xU 0 > yfto - yU o, 



(3.9) 


and comparison of (3.3), (3.9) then confirms our previous result (3.6). 

In practice the numerical solution of (2.9) and indeed the reduced large 7 problem, (3.2), is not straight- 
forward. The reason why there is a difficulty with the numerical solution of (3.2) was first discussed by 
Davey (1961) and later in more detail by Schofield and Davey (1967). In order to see what this difficulty is 
we consider the large ( limit of the equations to determine (u, fl, iD, f7) in (2.9) with N — 0. Suppose that 
for £ >> 1 we write 

(&,fl,iD, U) = (u + ,v + ,u> + + tUoo, 1+ U+) 


7 


where u + , v + , u>+ and 17+ are small and for simplicity we assume that E7 + is real. It is an easy matter to 
show that the linearized equations for (u + , v + , u> + , 17+) can be reduced to 


(!*-"”£) “ + = ^“ + . 


(3.10o,6) 


and U+, can be found in terms of u+,t> + . Thus for large £ we can write 

u+ = Aie™ 1 ^ + + A 3 e m *t } 

v + = B x e w ~<. 

Here m X} 7713 , m 3 are given by the values of m which satisfy 

2 m- Woo ± y/ w lo ± 87 , mr < 0 , 


(3.11) 


W- 


and we note that two of the required values will be complex when 7 > — ^ . Thus we have five independent 

o 

constants A Xi A 3 , A 3 and B 1 which can be iterated upon in order to satisfy the conditions at £ = 0 . 
It follows that there will be a continuum of solutions of (2.9) since there is no reason to ignore any of the 
decaying exponential solutions. Clearly it is to be expected that io*, will vary in this continuum so the 
solutions can be labeled by the size of the inflow towards the disc at infinity. This situation persists in the 
large 7 limit where 10 ^ ~ 7 ; Schofield and Davey (196^ argued that the solutions should in this case be fixed 
by discarding the slowest decaying exponential. Whilst it is certainly true that this fixes the solution, there 
is no basis for making such an assumption. Having made that assumption Schofield and Davey concluded 
that, in our notation, w = 87 . Interestingly enough we shall see later that this choice of w « fixes the 
boundary between linearly stable and unstable solutions of (2.9). We now present results obtained for the 
system (2.9) in the steady case N = 0 . 

As mentioned above, at large values of ( we have five constants, A*, B Xl A 2 , ^ 3 , 11/00 at our disposal 
once the constant 7 has been fixed. In our calculations we fixed 7,1000 and iterated on the remaining four 
constants until the required boundary conditions at ( = 0 were satisfied after integrating the differential 
equations in (2.9) from a suitably large value of ( to the origin. This integration was carried out using a 
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fourth order Runge-Kutta scheme or a compact finite difference scheme We restricted our attention to the 
cases 7 = .02, .1, .5 and in Table I we show the values of tt^O), ^(Oj^Ef^O) obtained for the different values 
of 7,iyoo shown. We see that for each of the values of 7 there are values of w ^ greater and less than 87. For 
each of the values of 7 used we were able to find solutions of (2.9) only for Wqo less than some critical value. 
Thus for example when 7 = .02 we were unable to find solutions of (2.9) for Wqq greater than— .355. Our 
calculations suggested that this minimum value decreases when 7 increases. In Figures (3. 1-3.4) we show the 
functions tl, 9, t2), U for different values of 7,100, Further solutions of (2.9) were obtained by fixing iWoo and 
varying 7. In particular we obtained solutions for the case w ^ = 87 and our results are shown in Figures 
(3.S-3.6). We recall that Davey obtained solutions of the large 7 problem in a different context and that his 
results were obtained by neglecting the slowest decaying exponential solution at large (. Having made that 
approximation Davey found numerically that w = 87 to the numerical accuracy of his calculations. Thus 
we expect that our results in Figures (3. 5-3. 6) should reduce to those of Davey at sufficiently large values 
of C. In fact we have in these figures shown Davey’s results expressed in our notation and we see that our 
results approach (3.7) for large 7. However we stress at this point that there is no reason why the solutions 
of (2.9) obtained by rejecting a particular decaying exponential solution of that system at large ( should 
be preferred, we hope to shed some light on the selection mechanism for the different solutions later in this 
paper. 

4 Time-dependent solutions of the interaction equations 

We shall now discuss the solutions we have obtained for the unsteady form of the interaction equations. 
We restrict our attention to the case when the stagnation point flow vanishes at infinity. This means that 
we are in effect discussing the finite amplitude instability of Karman’s solution. In the Appendix we shall 
give a limited discussion of the more general problem with 7^0. The appropriate simplified form of (2.10) 
is found by setting 7=1 and applying the conditions U = 0,£ = 0, 00; we obtain 
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U ^ — Ut = 2uTJ + 

Q« - Or = O 3 + |£7| 2 - 0 3 + iDGf , 
0 << - ©T = 200 4 - tf> 0 <, 


0= ® = U = 0,0 = 1,< = 0, 
0 = 0=17 = 0, £ = oo, 


(4.1) 


o = o(0,v = o(0,^ = ^(0,r = oJ 

The above system is parabolic in T and can be solved by marching forward in time from T = 0; we 
note here in passing that i D cannot be specified arbitrarily at T = 0 and must be deduced from u via the 
equation of continuity. For large values of T the solution of (4.1) will approach Karman’s solution if that 
flow is stable. We can therefore regard (4.1) as the nonlinear initial value instability problem for Karman’s 
rotating disc flow. However, we should bear in mind that (4.1) describes only finite amplitude disturbances 
with azimuthal wavenumbers ±2. 

In the first instance we restrict our attention to small initial perturbations from Karman’s solution, we 
therefore write 


(fi,€, U) = ( 0 , 15 , 0 ) + (u*,w*,£ 7 *) 


where (tt, 0) is Karman’s solution and u* etc are small. We now substitute 


(0, 0, tD, U) = (O + 0, 0 + v, ti> + w, U) 

in (4.1) and linearize to obtain the following decoupled problems: 

U i( - U T = 2617 + &U ( , ' 

17 = 0,£ = 0,oo, > 

U = U\T= 0 , 


(4.2) 
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and 


- Ut = 26u - 2Vv + G>6( + wfy, ' 
V(( — Vt — 26€ + 2fiu -f WV( + wi>(, 


2u+ W( = 0, 

u=C = iIi=0,C = 0, 
u = v = 0, C = oo, 

(«,«) = (u>*), T = 0. 


> 


(4.3) 


It should be pointed out that in the above equations we have in effect assumed that 0(17*) ~ 0(u*) ~ 
O(v’); a modified form of the equations can be derived when 0(U*) ~ O(u’)^ . In that case a nonlinear term 
|17| 3 must be inserted into the right hand side of the u equation. This particular case would be important 
only if the U problem were unstable. We can integrate (4. 2), (4.3) formally by taking a Laplace transform 
in time, when the transform is inverted the nature of the solution will depend crucially on whether either of 


the eigenvalue problems 

y" - ay = 2 Gy + Gn/ , 

y(o) -y{ oo) = o, 

0r y" - ay = 26j/ - 2 Gz + Gry 1 + *G', ^ 

z" -az = 26 z + 2yQ + rbz' + *«', 

► 

2y + x' = 0, 

*(0) = y(0) = z(0) = y(oo) = z(oo) = 0, , 

has an eigenvalue a with positive real part. Thus we shall now discuss these eigenvalue problems. 


(4.4) 


(4.5) 


The above eigenvalue problems were solved numerically; no unstable eigenvalues of either system were 


found so we conclude that the Kar man’s solution is stable to small amplitude perturbations of the type 


discussed here. In fact no discrete stable eigenvalues were obtained either, this is because both eigenvalue 
problems have a continuous spectrum over part of the plane ay < 0. The origin of this continuous spectrum 
can be seen from (4.4) by taking ( » 1. We see then that the two exponential solutions of the equation 
for y both decay if a is within the parabola a T = -a?/w ^ where Woo = G;(oo). Thus in this region we 
can always find a solution of (4.4) by combining the two independent solutions for y to satisfy the required 
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condition at the wall. A similar continuous spectrum can be seen to exist for the system (4.5); we expect 
the continuous spectra play an important role in the initial value problems (4.2),(4.3). Indeed since there is 
apparently no discrete spectrum associated with (4.4), (4.5) it is clear that the initial value problem must in 
some sense be described completely by the continuous spectrum. 

The initial value problem can be solved by taking Laplace transforms and inverting for particular forms 
of the initial perturbation, these inversions cannot in general be carried out analytically but their large time 
behaviour can be approximated asymptotically in a routine manner. Rather than use the Laplace transform 
method we Bhall instead look directly for the large time behaviour of (4.2), (4.3). We shall in fact restrict 
our attention to (4.2), but a similar approach can be used for (4.3). 

Suppose then that w Q 0 denotes the limiting value of fi) at large values of we choose to express ^ in 
the form 

& = (4. 6 ) 

where the function N then satisfies the modified equation 

+ (u>oo - ft )N( + (-^00 +™oo) jy = N t + 2fiJV. (4.7) 

At this stage we assume that all the exponential time dependence of the disturbance has been taken out 
by the substitution (4.6) so that N has only an algebraic dependence on T. It is well-known for the heat 
equation that the similarity variable ^ essentially replaces < when the the initial value problem is solved. 
Here the situation is slightly more complicated and we must seek a solution of (4.7) for £ = 0(1) and 
( = O(T^) for large values of T. Thus for ( = 0(1) we write 

N = TNo(C) + ---, 

where j is to be determined and JVq satisfies the ordinary differential equation 

K' + ( w oo - 0>)K + = 2 ujv 0 . (4.8) 

This equation must be solved subject to N 0 = 0, C = 0 and the resulting solution will then have N 0 ~ c( 
for large (, here c is an arbitrary constant which can be set equal to unity but whose actual value depends 
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on the form of the initial disturbance. Now let us find the required form for N in the upper region, before 
doing so we note that the constant j is left unknown at this stage since it plays no role in the zeroth order 
solution in the lower region. In the upper layer we write 

N = T #+ *Nb(x) + '”» ( 4 - 9 ) 

where x » the similarity variable Note here that the matching condition with the lower layer now 

rcquircB that for small x> ~ X* 

The equation to determine No is found to be 

tfo+!*o-(i + ^)tfo = 0. (4.10) 

In order that the disturbance decays to zero at large values of x we must insist that N 0 behaves like the 
exponentially decaying solution of the above equation; since it must also go to zero like x for small x we can 
show that the required solution is 

. _ 2*e^(-n-j,2-*x) (4. 11) 

V'(-n-lO) 

Here n is an odd integer, U(a, x) is the parabolic cylinder function whilst the constant j has been choosen 
to satisfy the matching condition at x = 0, this gives 



It follows that the solution will be dominated by the j = -§ eigensolution for large enough values of the 
time, the overall amplitude of this and the other decaying modes can only be determined by solving the 
initial value problem. In order to verify the above large time behaviour of the solution of the linearized 
perturbation equation for U we integrated (4.2) forward in time from T = 0 for the three cases: 

CASE a U* 

CASE b U* = < cos(e ~ ( , 

CASEc U* = 

We note that it is sufficient for ub to consider real initial data for if we are solving (4.2). The results 
we obtained are shown in Figure (4.1); in order to pick out the dominant exponential decay factor of U 
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we have plotted {log(T$ U' (0 ,T))t ■ On the basis of our discussion above we see that this quantity should 
2 

tend to — ~ — .2 for large T. We see that each of the above cases leads to results consistent with our 

predictions. A similar analysis to that given above for (4.2) can be given for (4.3), again the outcome is that 
a two-layer structure is required to describe the large time behaviour of the disturbance, furthermore the 
functions u, v are found to decay exponentially for large T with the same decay rate as that found above. 

We shall now report on some calculations carried out for the full nonlinear problem (4.1) with initial 
conditions 

*=0 = 0, U = (4.13) 

for 6 — .35, .45, .55, .65 . The results obtained for the initial conditions given above are typical of those we 
have found for a wide range of disturbances. The results we found are shown in Figure (4.2a,b,c,d,) where 
we have shown the growth rates (log /” (log /“ £T a dc) , (log 7 a d<) ,(log(U i (0,T).Ti) T 

. For the two smaller values of the amplitude constant 6 we see that the disturbance decays to zero so 
that Karman s solution is stable, note that Figure (4.2a) for 6 — .35 confirms the linear decay rate, -.2, 
at sufficiently large values of T when the disturbance field is now small. The calculations for the three 
larger values of 6 demonstrate that Karman *s solution is subcritically unstable. At a finite value of T our 
calculations encountered a singularity and could not be continued further. We did, of course, check that 
the singularity remains when the ( and T step lengths were decreased. In Figures (4.3a,b,c,d) we show the 
velocity field for the case T = 3.3 , 6 = .55. Further calculations for other initial disturbances confirmed the 
threshold type of response described above. 

Thus we have found that Karman ’s rotating disc solution is unstable to finite amplitude n = ±2 
modes wheras in the linear regime we have stability to this type of disturbance. We have made no attempt 
to investigate the dependence of the required threshold amplitude of the instability for different initial 
disturbance profiles. For the case discussed above instability occurs when the disturbance velocity field is of 
size comparable to the unperturbed state. Thus this particular type of disturbance is unlikely to be present in 
an experimental investigation so that it is unlikely that the this disturbance would cause transition. However 
it is not unreasonable to expect that a detailed investigation of a more general class of initial conditions 


14 



would isolate much more dangerous disturbances likely to cause instability in an experimental facility with 
moderate background disturbances. 


5. The structure of the Singular solutions 

Suppose that the singularity encountered in the numerical calculations occurs at time T = T. For 0(1) values 
of £ then we expect that the in viscid terms in (2.9) will dominate the flow; thus in the neighbourhood of T 
we expect that V ~ U ~ iD ~ 0 (T — We therefore expand the velocity components as 

VO 


U - 


9 = 
t& = 


T-T 

”»(C) 

t-t 

MC) 

T-T 


+ • - - | 
+ « •• > 


(5.1a, 5, c, d) 


u = Ml + ... . 

T-T 

We can then show from (2.9) that the system to determine (u,,u,,iu,, Z7,) is 

U, = u, 3 — v, 3 + 7 3 Uj + w,u,( , 
v f = 2u t v $ + , 

U t = 2U $ u $ + w t U t { , 

w t £ + 2u, = 0 

u>i - 0, C = 0 > 

v, = v # = U 9 = 0, ( = °o . 

Since the boundary conditions and differential equations for U, and v 9 are the same we can write 


(5.2a, 5, c,d) 


(5.3a, b) 


U a = Av $ 


where A is a constant. In fact the equation for v, is easily integrated to give 


v t — Qwgefoo^'k , 


where Q is a constant and w soo is ty,(oo). It is sufficient for us to choose Q = 1 and write the u, equation 
in the form 

u, - u, 2 + Sv 9 2 + w s v t( (5.4) 
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when 6 is a constant to be specified and it lies in the range (- 00 , oo). The case 6 = 0 corresponds to U 9 = v 9 
and in that case the equations for u Bt w 9 can be integrated directly to give 

u>t = U’.oot'U - l ) 3 , 

c 


u t — inu 9 = 1 — 


(5,5a, b) 


w 9 


(5.6a, b) 


It is easy to show that (5,5b) then defines a single- valued function u 9 = u 9 (£) satisfying 

u 9 ~ 1, ( — ♦ 0 

u t ~ , (-4 00 , 

For small values of £ the above solution must be matched onto a wall layer solution valid for £ = 0(T — T)“i . 
However we will now return to the more general case 6 ^ 0; it is easy to show that w t satisfies 


/2 


" 


w 9 ~ — ^ 6 


• OO -|- WgW t 


(5.7) 


and for large values of ( we then find that 


w, w too + Be-— , 


(5.8) 


where B is an arbitrary constant which is fixed by the condition u?,(0) = 0. For small values of ( we can 
show that Wg , v 9 have the forms 

(5.9a, 6) 

w, ~ 2 C + a£$ H 

where a, b are constants fixed by the choice of the arbitrary constant B appearing in the large £ form of 
Wg . In Figure (5.1) we show the function w 9 obtained by integrating (5.7) for three values of 6. In these 
calculations we have fixed w 9oo = — 1. It remains for us to show how the inviscid forms (5.9) for ( « 1 
can be matched onto a wall layer solution valid for ( = 0(T - ?)”$. If we define the similarity variable 

C 


V = 


(T-T) 


y then, in view of (5.9), we see that the appropriate expansions 


are 


An) M 

tl — — — h • • • , 

(T - T) ’ 

(t- fy 
w°(v) , 

w = i - i -r 4 — • 

(T-Ty 


(5.10a, b, c) 
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The equations to determine u°, v°, to 0 are found to be 

u°" + tj/2u 0 ’ + = u 0 ’ + w°u° , 

v°” + t}/2v 0 ' + ^t)° = w°v 0 ' + 2t»°ti 0 , 

w 0 ' + 2u° = 0, > ( 5 - U ) 

U° = V° = W° = 0 , 7} — 0 

U° ** 1, w° ~ 7] } 77 — ► OO, , 

From the asymptotic forms of the equations for u°, uj° at large 77 we find that there are three free constants 

which can be chosen so that the required conditions at 77 = 0 are satisfied. Thus we conclude that the viscous 
wall layer at ( = 0 (T - T)~b smoothes out the singular inviscid solutions so as to satisfy the boundary 
conditions. We conclude then by noting that the singularity structure given above is of course a singularity 
of the full Navier StokeB equations; in the physical situation we would anticipate that as the singularity 
develops the flow becomes unstable to disturbances with a more general 9 dependence. Presumably the 
growth of these other modes destroys the singularity structure and the Navier Stokes equations remain valid. 

6. Conclusion 

We have seen in the previous section that sufficiently large initial perturbations to Karman’s rotating 
disc flow lead to the development of a singularity of the Navier Stokes equations. We note that the structure 
of the singularity we have described does not depend at zeroth order on the nature of the underlying basic 
state. Therefore it is likely that the singularity structure we have found could be set up by sufficiently large 
amplitude perturbations to other basic states. 

If we assume that the structure we have found can be induced experimentally then a question of some 
importance is that of how the singularity can be controlled through the Navier Stokes equations once it has 
begun. Since no terms in the Navier Stokes equations have been neglected in the analysis leading to the 
singularity it might appear that the singularity must be controlled by an alternative set of field equations. 
However we do not believe that is the case, more precisely we believe that once the singularity has begun to 
develop the velocity profiles associated with it will be massively unstable to inviscid modes with azimuthal 
wavenumbers ^ ±2; these modes will then grow and prevent the further development of the singularity. 
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The latter conjecture could of course be verified by Navier Stokes simulations, we do not attempt such an 
investigation here. 

The typical size of initial amplitude associated with (4.13) required to cause the development of a singu- 
larity was found to be comparable with a typical basic state velocity. Thus experimentally the breakdown we 
have described could only be provoked by a disturbance amplitude unlikely to be present in an experimental 
situation. However it would be very surprising if the initial amplitude required to induce the singularity 
could not be significantly reduced by allowing a much more general initial perturbation. A possible way of 
isolating the most dangerous type of initial perturbation would be to formulate the Energy Stability problem 
associated with the interaction equations. 

Another open question following our discussion in Section 3 is that of which, if any, of the different 
equilibrium states for non zero 7 are most physically relevant. Our linear instability analysis in the Appendix 
strongly suggests that all those equilibrium states with w ^ < 87 are not relevant physically since they are 
linearly unstable. It remains an open question as to whether nonlinear effects are able to further reduce the 
class of physically realizable flows. 

One of the authors, PH, wishes to thank SERC and AFOSR for partial support of this work. 
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Appendix* The linear instability problem for the steady states with 7 ^ 0. 

Here wc shall discuss the instability of the steady equilibrium solutions of (2.9) with N = 0. We denote 
this steady state by (ft, 5), U) and suppose that we consider the instability of this flow to a time-dependent 
small perturbation (ft, 5, w } U) such that 

(ft, ft, IT) = («*,«*, T = °- ( A1 ) 


From (2.10) the linearized perturbation equations for (ft, ft, <b y U) are 

4 - 2 vU 4- 2 uO 4- 4- 

ur 4- 2tt+2'y*UU - 2ftft 4- ft>ft<; 4- ft)ft< = 

vt 4- 2ftft 4- 2ftft + ®ft{ 4- ft>ft( = ft«, 

2ft 4- = 0, 

ft=ft = ft) = 17 = 0, C = 0, 


(A2) 


subject to 


(A3) 

u = v = U = 0 , C = 00. 

and (Al). Again the initial value problem can be solved by a Laplace transform in T but the form of (A2) 
means that the resulting ordinary differential equations in ( must be solved numerically. Following our 
approach in §4 we consider the eigenvalue problem a = < 7 ( 7 ) obtained by replacing dr by a in (A2) and 
applying (A3). The structure of the disturbed flow at ( = 00 is then found by letting ( — » 00 in (A2) after 

replacing dr by <7 and <9< by m; the appropriate equations to determine m are 

(m 3 - Woom- cr)U = 2u, 

(m 3 - Woom - a)u = 2yU, (A4a, b, c) 

(m 3 - Woom - a)v = 0. 

Thus m is given by 

2m = u>oo ± V w lo + 

(A5a,6) 


2m = w 0O ± \fw\o + 4[tr ± 27 ]. 

<7? , , 

The first of these equations corresponds to (A4c) and leads to the continuous spectrum a r < — a found 

^OO 

in §4, However (A5a,b) lead to the continuous spectra, 


< r r < ±27 5-. 


( A6a, b ) 
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The positive root corresponds to the case when u = —7^ in (A4a,b) and, surprisingly, we obtain an unstable 
continuous spectrum for 7 > 0. However, it remains to be seen whether this unstable continuous spectrum 
can induce a physically relevant exponentially growing solution. In order to answer this question we shall 
now seek a large time solution of (A2); the structure we choose is based on the assumption that at large 
times the unstable spectrum (A6a) effectively dominates the flow. The first step in our solution procedure 
is the substitution 

= ( u +,v + ,«; + ) tr+)exp|^.C-^T+27rJ. (A 7 ) 

The functions u+,t; + etc. are then found in the regions where ( = 0(1), £ = O(T^). In the region where 
( = 0(1) we write 

(u+,v + ,ii/ + ,tf + ) = T j (u$v£,wZ } U£) + --- 

where uj etc. depend only on £. The problem for these functions is obtained from (A2) and solved subject 
to 

^0 = Vo = w o = u o = 0 C = 0 

and the unknown quantities uj (0), Vq ’(0), 10^ are chosen such that 

u o ~ 7 Cj &o ~ —() v o — 1 ► 0 exponentially £ — ♦ 00 

The constant j again remains unknown in the solution of the ( = 0(1) problem; however an investigation of 
the problem with £ = O(Ti) shows that matching with the above solution can only be achieved if j is again 
determined by (4.12). Thus we conclude that the steady state solutions with Wqq > 87 are linearly unstable 
and so not likely to be observed. 


20 



TABLE I 


^OO 

m 

*(0) 

7 ^( 0 ) 

- 0.35500000 

0.51152874 

- 0.60931513 

0.01281756 

- 0.39000000 

0.51089963 

- 0.61151223 

0.01015261 

- 0.41000000 

0.51080189 

- 0.61208108 

0.00945387 

- 0.45000000 

0.51069253 

- 0.61291648 

0.00842226 

- 0.50000000 

0.51062290 

- 0.61366951 

0.00748787 

- 0.55000000 

0.51058552 

- 0.61424223 

0.00677497 

- 0.60000000 

0.51056418 

- 0.61470076 

0.00620302 

- 0.65000000 

0.51055185 

- 0.61507974 

0.00572962 

- 0.70000000 

0.51054493 

- 0.61540001 

0.00532912 

- 0.75000000 

0.51054144 

- 0.61567522 

0.00498469 

- 0.80000000 

0.51054015 

- 0.61591483 

0.00468459 

- 0.85000000 

0.51054034 

- 0.61612571 

0.00442034 

- 0.90000000 

0.51054151 

- 0.61631296 

0.00418557 

- 0.95000000 

0.51054334 

- 0.61648051 

0.00397542 

- 1.00000000 

0.51054560 

- 0.61663143 

0.00378607 

- 1.05000000 

0.51054815 

- 0.61676816 

0.00361447 

- 1.10000000 

0.51055088 

- 0.61689265 

0.00345817 

- 1.15000000 

0.51055370 

- 0.61700654 

0.00331516 

- 1.20000000 

0.51055657 

- 0.61711115 

0.00318377 

- 1.25000000 

0.51055945 

- 0.61720760 

0.00306261 

- 1.30000000 

0.51056231 

- 0.61729682 

0.00295050 

- 1.35000000 

0.51056513 

- 0.61737962 

0.00284646 

- 1.40000000 

0.51056790 

- 0.61745666 

0.00274962 

- 1.50000000 

0.51057324 

- 0.61759580 

0.00257472 


- 1.60000000 

- 1.70000000 

- 1.80000000 

- 0.51000000 

- 0.55000000 

- 0.60000000 

- 0.65000000 

- 0.70000000 

- 0.75000000 

- 0.80000000 

- 0.85000000 

- 0.90000000 

- 0.95000000 

- 1.00000000 

- 1.10000000 

- 1.20000000 

- 1.30000000 

- 1.40000000 

- 1.50000000 

- 1.60000000 

- 1.70000000 

- 1.80000000 

- 1.90000000 

- 0.60000000 

- 0.80000000 

- 1.00000000 

- 1.50000000 


0.51057830 

0.51058306 

0.51058753 

0.52025270 

0.51846172 

0.51796137 

0.51779186 

0.51775162 

0.51777457 

0.51783140 

0.51790706 

0.51799309 

0.51808449 

0.51817818 

0.51836538 

0.51854617 

0.51871723 

0.51887748 

0.51902693 

0.51916605 

0.51929552 

0.51941606 

0.51952843 

0.68197809 

0.68831733 

0.69755716 

0.71653343 


- 0.61771805 

- 0.61782635 

- 0.61792299 

- 0.60663520 

- 0.61209594 

- 0.61511910 

- 0.61725674 

- 0.61895371 

- 0.62037081 

- 0.62158902 

- 0.62265639 

- 0.62360449 

- 0.62445542 

- 0.62522548 

- 0.62656993 

- 0.62770915 

- 0.62869011 

- 0.62954560 

- 0.63029943 

- 0.63096949 

- 0.63156952 

- 0.63211033 

- 0.63260051 

- 0.66231775 

- 0.71390532 

- 0.73310149 

- 0.75988876 


0.00242100 

0.00228479 

0.00216323 

0.03934270 

0.03382640 

0.03074344 

0.02852767 

0.02674543 

0.02524172 

0.02393845 

0.02278891 

0.02176219 

0.02083638 

0.01999520 

0.01851929 

0.01726178 

0.01617416 

0.01522217 

0.01438070 

0.01363075 

0.01295762 

0.01234970 

0.01179769 

0.20129254 

0.16251868 

0.14469129 

0.11597686 



.5 

- 2.00000000 

0.73020710 

- 0.77552448 

0.09744803 

.5 

- 2.50000000 

0.74055001 

- 0.78014433 

0.08420084 

.5 

- 3.00000000 

0.74852020 

- 0.79390800 

0.07420019 

.5 

- 3.50000000 

0.75489053 

- 0.79985908 

0.00030331 

,5 

- 4.00000000 

0.70008093 

- 0.80457555 

0.00003820 

.5 

- 4.50000000 

0.70441004 

- 0.80841092 

0.05482440 

.5 

- 5.00000000 

0.70800312 

- 0.81159300 

0.05045039 
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Figures 


Figure (3.1) The function tl(C) for different values of 7,14)00. 

Figure (3.2) The function c(f ) for different values of 7,14)00. 

Figure (3.3) The function Ti(f for different values of 7, u>oo • 

Figure (3.4) The function 717(f) for different values of7,uj 00 . 

Figure (3.5) The shear o'(0) for = 87. 

Figure (3.6) The shear yU'(0) for w = 87. 

Figure (4.1) The growth rate (logT$£7'(0, T))t for 
(a) IT * = (b) U* = Ccosfe-<, (c) U* = fe-«" 3 ) a . 

The predicted growth rate at large T is — ~ —.2. 

Figure (4.2a) The growth rate ai = (log /“ u 3 df ) T for 6 = .35, .45, .55, .65. 

Figure (4.2b) The growth rate cr a = (log / 0 °° v 2 dQ T for 6 = .35, .45, .55, .65. 

Figure (4.2c) The growth rate u 3 = (log /” U 2 d() T for 6 = .35, .45, .55, .65. 

Figure (4.2d) The growth rate <74 = (logT$ 17'(0, T)) t for 6 = .35, .45, .55, .65. 

Figure (4.3a) The velocity field fi(f ) for 6 = .55, T = 3.3. 

Figure (4.3b) The velocity field €(f) for 6 — .55, T = 3.3. 

Figure (4.3c) The velocity field tD(f) for 6 = .55, T = 3.3. 

Figure (4.3d) The velocity field 17(f) for 6 = .55, T = 3.3. 

Figure (5.1) The function u,(f) for lOoo = -1,7 = 0, .5, -.5. 
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Figure (3.5) The shear u’(0) f° r 



asympotic result 




Figure (3.6) The shear 7^'(°) for w 




Figure (4.1) The growth rate (log T*U'( 0, T))t for 

(a) U * = (e~<\ (b) U* = Ccos<e-<, (c) U* = (e~«- 

2 

The predicted growth rate at large T is — —.2. 
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for 8 = .35, .45, .55, .65. 
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Figure (5.1) The function uj#(C) for 
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